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Abstract 

We present some piloting calculations of the short-range correlation coefficients 
for the light and heavy baryons and masses of the doubly heavy baryons Hqq/ 
and Oqq/ {Q,Q' = c,b) in the framework of the simple approximation within the 
nonperturbative QCD approach. 
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1. Introduction 



The observation of meson by the CDF collaboration [[T| opens a new 
direction in the physics of hadrons containing two heavy quarks. Presently 
at the LHC, 5-factories, and the Tevatron with high luminosity, several 
experiments have been proposed, in which there is a possibility to identify 
and study hadrons containing two heavy quarks, like doubly-charm baryons 
{ccq) or baryons {hcq) with charm and beauty^. In the more distant future 
the next generation experiments with high bottom quark production rate 
will provide excellent possibilities for the study bottom baryons and their 
decays. In view of this project, it is important to have safe theoretical 
predictions for heavy baryon masses as a guide to the experimental search 
of these hadrons. 

A number of authors [0|-[lT^ have already considered baryons containing 



two heavy quarks in anticipation of future experiments which may discover 
these particles. In most of these works, however, theoretical predictions 
are somewhat biased by the introduction of the additional dynamical as- 
sumptions and supplementary dynamical parameters like constituent quark 
masses in addition to the only one parameter really pertinent to QCD - 
the overall scale of the theory Aqcd- 

The purpose of this paper is to calculate the masses of the heavy baryons 
in a simple approximation within the nonperturbative QCD, developed in 
|T^-[|rBj . This method has been already apphed to study baryon Regge tra- 



jectories [1T5| and, very recently, for computation of magnetic moments of 
light baryons [0 . The essential point of this paper is that it is very reason- 
able that the same method should also hold for hadrons containing heavy 
quarks. In this work we will concentrate on the masses of doubly heavy 
baryons. As in [jT?! we take as the universal parameter the QCD string 
tension cr, fixed in experiment by the meson and baryon Regge slopes. We 

^ Here, and throughout this paper, q denotes a hght quark u or d. 



2 



also include the perturbative Coulomb interaction with the frozen coupling 
as{l GeV) = 0.39. The basic feature of the considered approach is the dy- 
namical calculation of the quark constituent masses rrii in terms of the 
quark current masses mf^\ This is done using the einbein (auxiliary fields) 
formalism, which is proven to be rather accurate in various calculations for 
relativistic systems. The einbeins are treated as the variational parameters 
which are to be found form the condition of the minimum of baryon eigen 
energies [jT^ . 



2. Formalism 

The starting point of the approach is the Feynman-Schwinger representa- 
tion of the 3q Green's function, where the role of "time" parameter along 
a quark path is played by the Fock-Schwinger proper time. The final step 
is the derivation of the cm. Effective Hamiltonian (EH) containing the 
dynamical quark masses as parameters. For many details see the original 
papers []T^-lTBj. 



Consider a baryon consisting of three quarks with arbitrary masses m^, 
z = 1,2,3. In what follows we confine ourselves to consideration of the 
ground state baryons without radial and orbital excitations in which case 
tensor and spin-orbit forces do not contribute perturbatively. Then only 
the spin-spin interaction survives in the perturbative approximation. The 
EH has the following form: 

3 / (0)2 \ 



where mf^ are the current quark masses and rrii are the dynamical quark 
masses to be found from the minimum condition (see Eq. (|2|) below). Since 
rrii ^ mf^ for light quarks, but nii ~ mf'' for heavy quarks, each light 
quark contributes to the baryon mass an additional mass rrii/ 2 (not rrii 
as in the ordinary nonrelativistic quark model) , whereas each heavy quark 
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contributes an additional mass ~ mj. The dynamical quark masses are 
evaluated from the equations defining the stationary points of the baryon 
mass Mb as function of rrii 

^^ = 0. (2) 

orrii 

Let Ti be the quark coordinates. The kinetic momentum operator Hq in 
Eq. (|T]) acquires the familiar form 

^ _ _ 

^ 2midrl 2m2dr2 2m'sdrl' 

V is the sum of the perturbative Coulomb- like one gluon exchange potential 
and the string potential. The Coulomb-like potential is 

where the factor 2/3 is the value of the quadratic Casimir operator for the 
group S'[/c(3). The string potential has been calculated in []T^ as the static 
energy of the three heavy quarks 

Ktring(»^l, ■^2, ^3) = Cri^niin, (5) 

where i^min is the sum of the three distances |rj| from the string junction 
point, which for simplicity is chosen as coinciding with the centre-of-mass 
coordinate Rem- 

3. Hyper Radial Approximation 

We use the hyperspherical formalism approach (for detail see original pa- 
pers []T^). In the hyperradial approximation (HRA) corresponding to the 



truncation of the wave function ^^({^"i}) by the component with grand or- 
bital momentum K = the three-quark wave function depends only on 
the hyperradius = + A^, where p and A are the three-body Jacobi 
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variables^, and does not depend on angular variables. The confining po- 
tential (|5|) has a specific three-body character. However, this potential as 
well as the Coulomb potential in Eq. (g) is smooth in the sense that the 
HRA (where only the part of the potential which is invariant under ro- 
tation in the six-dimensional space spanned by the Jacobi coordinates is 
taken into account) is already an excellent approximation. The HRA ne- 
glects the mixed symmetry components of the three-quark wave function, 
which appear in the higher approximations of the hyperspherical formalism 
10]. Introducing the reduced function^ x(^) = R^^'^'ip{R) and averaging 



V = K:+Kitring ovcr the six-dimcnsional sphere one obtains the Schrodinger 
equation 



d'x{R) 



+ 2^ 



E - W{R) 
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x{R) = 0, (6) 



8fiR\ 

where fi is an arbitrary parameter with the dimension of mass which drops 
off in the final expressions. The last term in (0) represents the three-body 
centrifugal barrier and W{R) is the average of the three-quark potentials 
over the six-dimensional sphere: 

W{R) = {V) = -^ + bR, (7) 

with 

2as 16 32 ^ 

The mass depending constants aij and 7jj are defined by Eqs. ( \A.'4 ) and 
( |A.i;j| ) in the Appendix. 

It is convenient to introduce a new variable x = R^/Jl, to eliminate an 
artificial dependence of Eq. (g) on /i, then the equation (§) becomes 

+ 2(e- U{x) - ^) x{x) = , (9) 



^for their definition see Appendix 

■^In what foUows we omit the value of if = to avoid subscripts. Note that the radiaUy symmetric 
component with if = is the dominant one in the three-quark wave function. 
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where 



U{x) = ^ + X. 10 



Since a ~ l/v^5 ^ (^^^ (^3)' eigenvalue ^ in 

does not depend on fi. 

4. The quark dynamical masses 

Equation (0) applied to the nucleon {m^q'' ~ 0) yields the dynamical mass 
ijiq of the light quark, and applied to the strange hyperons gives the 
strange quark mass nis. In the same manner application of this equa- 
tion to the charm and beauty baryons yields the constituent masses of c- 
and 6-quarks. In our calculations we use the same parameters as in [|22 



namely a = 0.17 GeV, a, = 0.4, mf^ = 0.009 GeV, m?^ = 0.17 GeV, 
m'^c^ = 1.4 GeV, and mf^ = 4.8 GeV. 

We solve Eq. using both the quasiclassical and variational solutions. 

The first approach is based on the well-known fact that interplay between 

the centrifugal term and the confining potential produces a minimum of 

the effective potential specific for the three-body problem. The numerical 

solution of (0) for the ground state eigen energy may be reproduced on 

a per cent level of accuracy by using the parabolic approximation for the 

effective potential [pO], [j2Tl. This approximation provides an analytical 

15 

expression for the eigen energy. The potential U {x) = U {x) + ^— ^ has the 
minimum at a point x = xq, which is defined by the condition U'{xq) = 0, 
i.e.: 

b .3 



XQ^{a^)xo- 15/4 = 0. (11) 



Expanding U {x) in the vicinity of x = xq one obtains: 

U{x) ^ U{xo) + ^U'\xo){x - xof, 
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i.e. the potential of the harmonic oscillator with the frequency 



u = yU"{xo). Therefore the energy eigenvalue is 

Eo^U{xo) + ^u;. (12) 

In Table 1 we show the dynamical masses rrii and the ground state eigenval- 
ues Eq for various baryons calculated using the procedure described above. 
Our values of light quark mass niq qualitatively agree with the results of 
|22|] obtained from the analysis of the heavy-light ground meson states, 



but ~ 60 MeV higher than those of [1T5|, []r7|. This difference is due to 
the different treatment of the Coulomb and spin-spin interactions. In [|T5 



both interactions have not been included and the light quark mass has been 
calculated from the fit of the mass of A (1232) where the Coulomb-like po- 
tential and the spin-spin interaction seem to balance each other. In [|T7[ 
the smeared spin-spin interaction for the light quarks has been included 
into Eq. (0) defining the dynamical mass of the light quark. In our calcula- 
tion as in []22| we include the Coulomb-like term, but neglect the spin-spin 
interaction. 

There is no good theoretical reason why dynamical quark masses need to 
be the same in different mesons and baryons. From the results of Table 
1 we conclude that the dynamical masses of the light quarks (w, or s) 
are increased by ~ 100 — 150 MeV when going from the light to heavy 
baryons. For the heavy quarks (c and b) the variation in the values of their 
dynamical masses is marginal. In Table 2 we compare the quark masses in 
Aq and Sg baryons with those calculated in in D and B mesons. One 
observes that the masses of the light quarks in baryons are slightly smaller 
than those in the mesons. The small variations in the values of rric and 
are within the accuracy of our calculations. 
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5. Correlation functions for the baryons 

For many applications the quantities {^JJ\6'^^\rj — ri)\ip) are needed. To 
estimate effects related to the baryon wave function we solve Eq. (0) by 
the variational method. We introduce a simple variational ansatz for x{^) 



X{x) = 2V2p^x'/^e-P"''" , (13) 

where p is the variational parameter, and the numerical factor is chosen 
so that J 'x^{x)dx = 1. The trial three-quark Hamiltonian admits explicit 
solutions for the energy, the wave function, and the density matrix: 



EQ^minE{p), (14) 
p 



where 



E{p) = {x\H\x) = 3/ - {a^) ■\\ll-P+ (Vv^) • ^^f • P''- (15) 

The density matrix (the correlation function) fijk{fij) in a baryon {ijk} is 
defined as: 

f^^k{r^J) = a% J \i^{a,jr,j, \j)fd^X,j (16) 



so that 



J f^Jk{r^J)d^r,J = JJ |^(p,^., \j)\^d^\jd^p^^ = 1 . (17) 

For the trial function ( [Tp fijk{fij) are evaluated expHcitly: 

/c \3/2 

si 



kk(r^i) = (^^j e-&l-«l , (18) 
with 

6j = 2^0 • (19) 

where /i^j is the reduced mass of the quarks i and j, and pq is to be found 
from the condition 

=0. 

P=Po 
8 



dp 



The expectation values fijk{fij) depend on the third or 'spectator' quark 
through the three-quark wave function. 
Let us define the quantities 



3/2 

) 



R^Jk = hk{0) = ( % ) (20) 



The corresponding quantity for a meson is denoted as Rij. The results of 
the variational calculations are given in Table 3 where for each baryon we 
show the variational parameters ^o, the quantities Rijk (in units of GeV^), 
and the average distances fij = ^ (r|^-) (in units of fm). The variational 
estimations of £"0 and quark dynamical masses do not differ from those 
shown in Table 1. 

Comparing the results of Table 3 with those of we obtain (see Table 



Rijk < -jRiji (21) 

and 

Rijk > Rijh if ruk < mi (22) 

Note, however, that if j are the light quarks, and the quarks k and / are 
the heavy, then Rijk ~ Riji {e.g., Rqqc ~ Rqqb) in agreement with the limit 
of the heavy quark effective theory. 

Our estimations for the ratios Rijk/ Rij agree with the results obtained 
using the nonrelativistic quark model or the bag model [j2^-[^ or QCD 



sum rules which are typically in the range 0.1 — 0.5. On the other 



hand, our result for disagrees with the one by Rosner |2B1 who es- 
timated the heavy-light diquark density at zero separation in A5 from 
the ratio of hyperfine splittings between E^, and baryons and B and 

Inequalities (^ij) and (|2^) were first suggested in from the observed mass splitting in mesons 
and baryons. 
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5* mesons and found Rqhu/Rid ^ 0-9 =t O-l? if the baryon splitting is 
taken to be m|* — ~ m|* — = (0.384 ± 0.035) GeV^, or even 
to Ruhd/ Rhd ^ 1.8 ± 0.5, if the surprisingly small and not confirmed yet 
DELPHI result ms* - m^, = (56 ± 16) MeV is used. 

From the results of Table 3 it follows that the correlation between two 
quarks depends on the third one. Note also that the wave function calcu- 
lated in HRA shows the marginal diquark clustering in the doubly heavy 
baryons. This is principally kinematic effect related to the fact that in 
the HRA the difference between the various fij in a baryon is due to the 
factor ^Jxf^j which varies between for rrii = rrij and ^Jl/rrii for 

mi <C rrij. 

In Table 5 we compare the short-range correlation coefficients in the doubly 
heavy baryons with those calculated in [j7[] using the pair- wise quark inter- 
action with the power-law potential, and in [jSj using the non-relativistic 
model with the Buchmiiller-Tye potential. 



6. Masses of doubly heavy baryons 

To calculate hadron masses we, as in [jT5|, first renormalise the string po- 



tential 

String ^ String ~l~ ^ ^ C'j, (^3) 

i 

where the constants Ci take into account the residual self-energy (RSE) 
of quarks. In principle, these constants can be expressed in terms of the 
two scalar functions entering covariant expansion of the bilocal cumulants 
of gluonic fields in the QCD vacuum [[T^, [T^. In the present work we 



treat them phenomenologically. To find Ci in (^3]) we assume, first, that 
the spin splittings of hadrons with a given quark content arise from the 
colour-magnetic interaction in QCD. Indeed, for the ground state hadrons 
the hadron wave functions have no orbital angular momentum, so tensor 
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and spin-orbit forces do not contribute. The second assumption is that the 
colour-magnetic interaction can be treated perturbatively []2^, ^ : 



^^spin = — t: — / , Rijk- (24) 

9 ^-^ mirrij 

Kj -J 

Because the colour-magnetic interaction between two quarks goes inversely 
as the product of their masses, the perturbative approximation improves 
as the quark mass increases. However, this approximation may not be 
good for the baryons containing light quarks^. In what follows we adjust 
the RSE constants Cj to reproduce the centre-of-gravity for baryons with 
a given flavour. To this end we consider the spin-averaged masses, such as: 

^^^±^ = 1.085 GeV, »d Ma + + 2M.- ^ 

2 4 ' V 7 

and analogous combinations for qqc and qqb states. Then we obtain 

C^ = 0.34GeV, = 0.19 GeV, ~ C5 ~ 0. (26) 

We keep these parameters fixed to calculate the masses given in Table 6, 
namely the spin-averaged masses (computed without the spin-spin term) 
of the lowest doubly heavy baryons. Our results are very similar to those 
obtained in [)7| using the pair-wise power-law potential. 



7. Conclusions 

In this paper we employ the general formalism for the baryons, which is 
based on nonperturbative QCD and where the only inputs are the string 
tension cr, the strong coupling constant Ofg, and two additive constants, 
Cq and Cs, the residual self-energies of the hght quarks. We present some 
piloting calculations of the dynamical quark masses for various baryons 
(see Table 1). The latters are computed solely in terms of a and ag and 
depend on a baryon. 



'"^Note that l/rmmj dependence in Eq. (24), if treated literally in the EH method, results in a collapse 
both in the pseudoscalar qq channel and the proton. That may be a signal of the Nambu-Goldstone 
phenomenon. 



11 



The second important point of our investigation is the calculation of the 
correlation functions for baryons. They are given, among the other things, 
in Table 3. We have also performed the calculations of the spin-averaged 
masses of baryons with two heavy quarks. One can see from Table 6 that 
our predictions are especially close to those obtained in []7[] using a variant 
of the power-law potential adjusted to fit ground state baryons. 
Evaluation of the spin-spin interactions requires inclusion of the K = 2 
hyperspherical components and/or more sophisticated treatment of the 
colour-magnetic interaction. We shall consider these calculations in the 
next publication. 
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Appendix 

Consider three quarks with arbitrary masses m^, ? = 1, 2, 3, and coordinates 
Vi. The problem is conveniently treated using Jacobi coordinates Pj ■ and 





rrii + rrij 




(A.l) 



where 




(A.2) 



Here fiij and Hij^k are the reduced masses 



mimj 





rrii + rrij 



nii + rrij + rrik 
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Altogether with the centre-of-mass coordinate Rem Jacobi coordinates de- 
termine completely the position of the system. The Jacobian of the trans- 
formation for the differential volume elements is 1, i.e., 

d^p^d^Xu = d^p^2d^X^2 = d^Pizd^Xii^. (A.4) 

The inverse transformations for the relative coordinates r^j = — Vj and 
'^k ~ R-cm 8.re 



1 / uimi + rrij) ^ r x ^\ 

r^J = — P^J, rk - Rem = - \ 7 7 A,^- . A.5 

aij y mk{mi + m2 + ma) 

The hyperradius is defined as = p'^j + Xfj and does not depend on 
the order of the quark numbering: 

R' = p\2 + A?2 = P32 + = P?3 + A?3- (A.6) 
Written in terms of Vij Eq. ( A-^ ) reads: 

R' = y\ r% . {A.7) 

In the centre-of-mass frame Rem = the invariant kinetic energy operator 
is written in terms the Jacobi coordinates ( |A.1| ) as 

^'-~2il[W'^dX')~~^[dR^^RdR^ ) ' ^^-^^ 

where K'^{Q) is angular momentum operator whose eigen functions (the 
hyperspherical harmonics) are 

K\n)Y^K] = -K{K + 4) V], (A.9) 

with K being the grand orbital momentum. In terms of Y|x] the wave 
function i/j{p, X) can be written in a symbolical shorthand as 



i;{p,X) = J2MR)Y[K]m. 



K 
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In the HRA K = and i/j = i/j{R). Note that the centrifugal potential in 
the Schrodinger equation for the radial function ipxiR) with a given K 

(X + 2)2- 1/4 

is not zero even for K = 0. For the reduced function x(^) = R^/'^iIj{R) 
one obtains after averaging the interaction over the six-dimensional sphere 
Eq. (D with 

+ ^string) (A. 10) 

One can easily see that the definition of A)) does not depend on the 

order of the quark numeration. 

In terms of the Jacobi coordinates the Coulomb and string potentials read: 



^string = Cr^Tijl'^ijI? (A. 12) 

where 



mk[mi + 777,2 + rn^) 



Using the relations 



valid for any pair (ij), one obtains Eqs. (0). 



14 



References 

[1] F.Abe et a/., CDF Collaboration, Phys. Rev. Lett. 81, 2432 (1998). 

[2] S.Capstick and N.Isgur, Phys. Rev. D 34, 2809 (1986). 

[3] M.J. Savage and M.Wise, Phys. Lett. B 248, 177 (1990). 

[4] M.J.White and M.J. Savage, Phys. Lett. B 271, 410 (1991). 

[5] E.Bagan, M.Chabab, and S.Narison, Phys. Lett. B 306, 350 (1993). 

[6] A.F.Falk, M.Luke, and M.J. Savage, Phys. Rev. D 49, 555 (1994). 

[7] E.Bagan, H.G.Dosch, P.Godzinsky, S.Narison, and J.-M. Richard, 
Z.Phys. C 64, 57 (1994). 

[8] D.Ebert, R.N.Faustov, V.O. Calkin, A.P.Martynenko, and V.A.Saleev, 
Z. Phys. C 76, 111 (1997). 

[9] S.S.Cershtein, V.V.Kiselev, A.K.Likhoded, and A.LOnishchenko, 
Mod. Phys. Lett. A 14, 135 (1999); Heavy Ion Phys. 9, 133 (1999). 

[10] R.Roncaglia, D.B.Lichtenberg, and E.Predazzi, Phys. Rev. D 52, 1248 
(1995). 

[11] J.C.Korner, M.Kramer, and D.Pirjol, Prog. Part. Nucl. Phys. 33, 787 
(1994). 

[12] A.K.Likhoded and A.LOnishchenko, |hep-ph/99 1^425 . 

[13] Yu.A.Simonov, Nucl. Phys. B 307, 512 (1988); Phys. Lett. B 226, 151 
(1989); 228, 413 (1989). 

[14] Yu.A.Simonov, Nucl. Phys. B 324, 67 (1989). 

[15] M.Fabre de la Ripelle and Yu.A.Simonov, Ann. Phys. (N.Y.) 212, 235 
(1991). 



15 



[16] A.Yu.Dubin, A.B.Kaidalov, and Yu.A.Simonov, Phys. Lett. B 323, 41 
(1994). 



[17] B.O.Kerbikov and Yu.A.Simonov, [hep-ph/ 000 1243 



[18] Yu.S.Kalashnikova, A.V.Nefediev, and Yu.A.Simonov, [hep- 
I ph/0103274| . 



[19] Yu.A.Simonov, Yad. Fiz. 3, 630 (1966); A.M.Badalyan and 
Yu.A.Simonov, Yad. Fiz. 3, 1032 (1966); F.Calogero and 
Yu.A.Simonov, Phys. Rev. 183, 869 (1968). 

[20] Yu.S.Kalashnikova, I.M.Narodetskii, and Yu.A.Simonov, Yad. Fiz. 46, 
1181 (1987). 

[21] B.O.Kerbikov, M.I.Polikarpov, and L.S.Shevchenko, Nucl. Phys. B 
331, 19 (1990). 

[22] Yu.S.Kalashnikova and A.Nefediev, Phys. Lett. B 492, 91 (2000). 

[23] D.B.Lichtenberg, Phys. Rev. D 54, 5887 (1996). 

[24] M.B.Voloshin and M.A.Shifman, Yad. Fiz. 41, 187 (1985) [Sov. J. 
Nucl. Phys. 41, 120 (1985)]; Zh. Exp. Teor. Fiz. 91, 1180 (1986) [Sov. 
Phys. JETP 64, 698 (1986)]. 

[25] B.Guberina, S.Nussinov, R.Peccei, and R.Riickl, Phys. Lett. B 89, 111 
(1979); N.Bilic, B.Guberina, and J.Trampetic, Nucl. Phys. B 248, 261 
(1984); B.Guberina, R.Riickl, and J.Trampetic, Z. Phys. C 33, 297 
(1986). 

[26] LBigi et a/., in B-decays, ed. S.Stone, 2nd edition (World Scientific, 
Singapore, 1994) p. 134; LBigi, |hep-ph/9508TO . 



[27] P.Conangelo and F.De Fazio, Phys. Lett. B 387, 371 (1996). 
[28] J.L.Rosner, Phys. Lett. B 379, 267 (1996). 



16 



[29] A.De Rujula, H.Georgy, and S.L.Glashow, Phys. Rev. D 12, 147 
(1975). 

[30] W.Kwong, J.L.Rosner, and C.Quigg, Ann. Rev. Nucl. Part. Sci. 37, 
325 (1987). 



17 



Table 1. The constituent quark masses uii and the ground state eigen 
energies Eq (in units of GeV) for the various baryon states. (The results 
obtained from the quasiclassical solution and from the variational one prac- 
tically coincide.) 



baryon 


mi 


m2 


ms 


Eq 


{qqq) 


0.446 


0.446 


0.446 


1.438 


(qqs) 


0.451 


0.451 


0.485 


1.414 


{qss) 


0.457 


0.490 


0.490 


1.392 


(sss) 


0.495 


0.495 


0.495 


1.370 


{qqc) 


0.519 


0.519 


1.502 


1.176 


(qsc) 


0.522 


0.555 


1.505 


1.157 


(ssc) 


0.589 


0.589 


1.507 


1.138 


{qqb) 


0.564 


0.564 


4.836 


1.057 


(qsb) 


0.567 


0.601 


4.837 


1.038 


{ssb) 


0.604 


0.604 


4.838 


1.019 


{qcc) 


0.569 


1.555 


1.555 


0.926 


(sec) 


0.604 


1.557 


1.557 


0.908 


{qcb) 


0.606 


1.616 


4.866 


0.783 


{scb) 


0.642 


1.618 


4.867 


0.765 


(qbb) 


0.636 


4.931 


4.931 


0.582 


{sbb) 


0.673 


4.931 


4.931 


0.565 
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Table 2. The dynamical quark masses for the ground state (qc), (sc), (g6), 
(s6) mesons and for the corresponding ground state baryons. 



State 


TUq 


rris 


rric 


rub 


{qc) 


0.529 




1.497 




(sc) 




0.569 


1.501 




{qqc) 


0.519 




1.502 




{qsc) 


0.522 


0.555 


1.505 






0.619 






4.84 


{sb) 




0.658 




4.842 


{qqb) 


0.564 






4.836 


{qsb) 


0.567 


0.601 




4.838 
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Table 3. Rijk in units of GeV^ and fij = y {r'^j) in units of fm. (The 
results are obtained from the trial functions ([T3D with the variational pa- 
rameters Po given in units of GeV^^/^^ in the first column. The results for 
light baryons are presented for completeness.) 



baryon 


Po 


^123 


^231 


^312 


fl2 


^23 


^31 


(qqq) 


0.472 


0.00564 


0.00564 


0.00564 


0.777 


0.777 


0.777 


(qqs) 


0.470 


0.00567 


0.00598 


0.00598 


0.775 


0.762 


0.762 


{qss) 


0.469 


0.00600 


0.00633 


0.00600 


0.760 


0.747 


0.760 


(sss) 


0.467 


0.00636 


0.00636 


0.00636 


0.746 


0.746 


0.746 


(qqc) 


0.454 


0.00626 


0.0113 


0.0113 


0.750 


0.615 


0.615 


(qsc) 


0.452 


0.00656 


0.0121 


0.0113 


0.738 


0.601 


0.615 


(ssc) 


0.451 


0.00688 


0.0121 


0.0121 


0.727 


0.602 


0.602 


(qqb) 


0.447 


0.00681 


0.0163 


0.0163 


0.729 


0.545 


0.545 


{qsb) 


0.446 


0.00711 


0.0176 


0.0163 


0.719 


0.531 


0.545 


{ssb) 


0.445 


0.00742 


0.0176 


0.0176 


0.708 


0.531 


0.531 


{qcc) 


0.439 


0.0116 


0.0296 


0.0116 


0.611 


0.447 


0.611 


(sec) 


0.438 


0.0123 


0.0294 


0.0123 


0.599 


0.448 


0.599 


{qcb) 


0.436 


0.0123 


0.0562 


0.0166 


0.599 


0.361 


0.541 


{scb) 


0.435 


0.0130 


0.0559 


0.0178 


0.587 


0.361 


0.529 


{qbb) 


0.438 


0.0181 


0.165 


0.0181 


0.527 


0.252 


0.527 


{sbb) 


0.437 


0.0194 


0.165 


0.0194 


0.515 


0.252 


0.515 



Table 4. The ratios of the squares of the wave functions determining 
the probability to find a light quark at the location of the heavy quark 
inside the heavy baryon and the corresponding meson. (The meson wave 
functions are taken from [[22].) 



Rued/ Rue 


Rscu/ Rse 


Rubd/ Rbd 


Rsbu/ Rsb 


0.436 


0.405 


0.373 


0.340 
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Table 5. Short-range correlation coefhcients Rijk- In the parentheses are 
shown the corresponding quantities calculated using the power-law poten- 
tial I?! . In the square brackets are shown correlation coefficients calculated 
using non-relativistic model with Buchmiiller-Tye potential. 



State 


^123 


^231 


^312 


{ccq) 


0.030 (0.039) [0.022] 


0.012 (0.009) 


0.012 (0.009) 


(ccs) 


0.030 (0.042) [0.022] 


0.012 (0.019) 


0.012 (0.019) 


{bbq) 


0.165 (0.152) [0.144] 


0.018 (0.012) 


0.018 (0.012) 


(66s) 


0.165 (0.162) [0.144] 


0.019 (0.028) 


0.019 (0.028) 


{bcq) 


0.056 (0.065) [0.042] 


0.012 (0.010) 


0.017 (0.011) 


(6cs) 


0.056 (0.071) [0.042] 


0.013 (0.021) 


0.018 (0.025) 



Table 6. Masses of baryons containing two heavy quarks 



State 


present 
work 
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E{qcc} 
^}{scc} 


3.69 
3.86 


3.70 
3.80 


3.71 
3.76 


3.66 
3.74 


3.61 
3.71 


3.48 
3.58 


E{qcb} 
n{scb} 


6.96 
7.13 


6.99 
7.07 


6.95 
7.05 


7.04 
7.09 




6.82 
6.92 


E{qbb} 
n{sbb} 


10.16 
10.34 


10.24 
10.30 


10.23 
10.32 


10.24 
10.37 




10.09 
10.19 



(a) ^Y\e additive nonrelativistic quark model with the power-law potential. 

Relativistic quasipotential quark model. 
'^^^ The Feynman-Hellmann theorem. 
(d) Approximation of doubly heavy diquark. 
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